A description of the yrast states in 24 Mg by the self-consistent 3D-cranking model 
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With the self-consistent 3D-Cranking model, the ground-state rotational band in 24 Mg is analysed. 
A role of triaxial deformation is discussed, in particular, in a description of the observed two I 7 * = 8 + 
states. 
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INTRODUCTION 

ifMg has been studied very well as a typical case of 
well-deformed light-mass nuclear systems. After self- 
consistent microscopic calculations of the non-relativistic 
and relativistic methods Q in the late 1980s, the nu- 
cleus is believed to have an axially-symmetric and pro- 
late shape in the ground state. Assuming the core of 
16 0, the ground state configuration is supposed to have 
eight valence particles occupying the d 5 / 2 orbits (four 
neutrons and four protons) . Since the Fermi level of this 
nucleus is in the beginning of the sd-shell (corresponding 
to the N = 2 harmonic oscillator shell), there are many 
open valence orbitals above the Fermi level, which may 
induce deformation, and consequently a collective rota- 
tion. Such a rotational band has been already identified 
in experiments (for instance, see Ref. 0,0)- A partic- 
ular interest is an existence of the two 1^ = 8 + states 
observed in experiment. These two states are energeti- 
cally close to each other (the difference is about 2 MeV) . 
After the study of Sheline et al. @], it is believed that the 
second 8 + belongs to the ground-state rotational band (g- 
band) . Valor, et al. analysed the g-band with the cranked 
Skyrme HF + BCS approach as well as the configura- 
tion mixing approach based on the generator coordinate 
method (GCM) However, their calculations assume 
axial symmetry for descriptions of intrinsic states. Up to 
7 7r = 4 + , they were able to reproduce the experimental 
data very well. The cranked mean-field calculation gave 
a fairly good agreement for 7 7r = 6 + , while the config- 
uration mixing approaches returned larger energies for 
the state. This is simply because these states at high 
spin were projected out from the non-cranked mean-field 
state. Interestingly, the cranked mean-field calculation 
for I 7 * = 8 + matches the observed energy of the first 8 + 
state, but the authors dismissed the agreement on the 
basis of the Sheline's analysis ||, and they speculated 
that the disagreement might mainly come from the tri- 
axial effects at high spin induced by the disappearance 
of the pairing correlation in their calculation. 



3D-CRANKED HFB METHOD 

Inspired by the study by Valor et al. Q , we performed 
the self-consistent cranking calculation allowing triaxial 
deformation in a self-consistent manner. As a new aspect 



in our study, not only ID-cranking but also 3D-cranking 
calculations were carried out. An advantage of the 3D- 
cranking model is that low- and high-K intrinsic struc- 
tures can be systematically studied 

The Hamiltonian used in our study reads 



H = Hn 



V, 



(1) 



where the first term describes the one-body part, which 
is the spherical Nilsson Hamiltonian in this study, while 
the second part is the two-body interactions, which 
is the pairing-plus-quadrupole force (so-called P+QQ 
force). The model space (valence space) to diagonalise 
the two-body part contains two major shells (N — 2, 3) 
in the spherical Nilsson model, in accordance with the 
Kummar-Baranger criteria for the P+QQ force §]■ The 
variational state is the Hartree-Fock-Bogoliubov ansatz, 
which is a generalised product state. With quasi-particle 
annihilation operators (3 q , the ansatz is expressed as 



|HFB> = n&|0), 



(2) 



where |0) is the vacuum for the canonical basis a m and 
a^. (In our case, the canonical basis correspond to the 
spherical Nilsson basis.) The canonical basis and the 
quasi-particle basis are connected by a unitary transfor- 
mation called the general Bogoliubov transformation [j| . 
The variational equation is derived for 

3 

5(RFB\H - ("iJi + ViBi) - ArAMHFB} = 0. 

2—1 r— p,n 

(3) 

In the above equation, Ji is the i-th component of the an- 
gular momentum operator (the index i takes i = 1,2,3), 
and N T describes the number operator for proton (r = p) 
and neutron (r = n) . Bi is an off-diagonal component of 
the quadrupole operator, defined as 



Bi 



2 (Of* 



Qkj 



(4) 



where the indices (i,j,k) should be placed in a 
cyclic manner. Each term with the Lagrange mul- 
tipliers (uji , /itj and A T ) is necessary to put con- 
straints in intrinsic states: (( Ji), (J2), (J3)) = 
(Jcos#, Jsinflsin^, Jsin0cos</>); (N T ) — N T ;{Bi) = 0. 
The last constraint is necessary so as to keep the orienta- 
tion of the angular momentum vector against the intrin- 
sic coordinate axes . The variational equation is solved 
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J 2 4 6 8 9 
P 0.37 0.38 0.37 0.35 0.31 0.29 

TABLE I: Evolution of the deformation in (3 as a function 
of the total spin J, which are obtained from the ID-cranking 
calculation. 

by means of the method of steepest descent. Details of 
the method are presented in Ref . . A deformed Nilsson 
+ BCS state is used for an initial trial state at J = 0. De- 
formation parameters and gap energies for the trial state 
are determined by referring to the calculations of the liq- 
uid drop model by the Moller and Nix • I n the present 
study, the deformation parameters for the ground state 
are chosen to be (/3, 7) = (0.347,0.0°) and the pairing 
gap energies are (A p , A n ) = (1.840 MeV, 1.962 MeV). 
All the physical quantities, such as energy, quadrupole 
moments (deformation), single-particle spin components, 
and gap energies, are self-consistently calculated under 
the above constraints in this framework. 

RESULTS AND DISCUSSIONS 

First of all, let us report the results from the 1D- 
cranking calculation. Despite a use of a simple separable 
interaction, the ground-state rotational spectra is repro- 
duced reasonably well (FigJIJ. As Valor, et al. com- 
mented in Ref.[6j, the gap energies disappear both for 
protons and neutrons before J = 6h (FigQ}. Triaxial de- 
formation gradually decreases from 7 = 0°. However, in 
J < Ah, triaxial deformation can be still regarded negli- 
gible. In other words, axial symmetry is kept fairly well. 
On the other hand, at high spin (J = 8,10ft), substan- 
tial triaxial deformation is formed (7 > —10°), and axial 
symmetry is clearly broken. It should be noted here that 
the convention for the quadrupole deformation parame- 
ters (/3, 7) in this study follows the Hill- Wheeler coordi- 
nates, which gives the opposite sign in 7 to the so-called 
Lund convention. A fact that the gamma deformation be- 
comes negatively larger implies that the nucleus is reach- 
ing the non-collective rotational state (7 = —60°), where 
the rotational axis corresponds to the shortest principal 
axis of the deformation. 

From this result, it is understandable why the crank- 
ing calculation by Valor, et al. was successful up to 
J — 6h and why the deviations from the experimental 
values become larger at higher spins. As mentioned ear- 
lier, they solved the HFB equation within the axial sym- 
metry constraint. 

In the HFB theory, the total spin is expressed as the 
sum of single-particle spins, that is, 

(Ji) = ^(j^) = ^(rimnPnm, (5) 
m mn 

where p is the density matrix and ji is the single-particle 
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FIG. 1: (top) Rotational energy obtained in the ID cranking 
calculation (solid line), and the experimental data (cross), 
(middle) Gap energies for protons and neutrons, obtained in 
the ID cranking calculation, (bottom) Gamma deformation 
obtained in the ID cranking calculation. 

angular momentum operator. The indices m and n are 
for the canonical basis, while the index i for the coor- 
dinate axes, that is, i — 1,2,3. Using this information 
of angular momentum, we can discuss nuclear structure 
with single-particle spins. In the Table [H] are the cal- 
culated main components of single-particle spins for the 
different total spin J. The result reflects a fact that 24 Mg 
is a N — Z nucleus, that is, the way of single-particle ex- 
citations is the same both for protons and neutrons. For 
low-spin members in the rotational band, the total spin 
consists mainly of the ds/2- The higher the total spin, the 
more the orbit is occupied. Therefore, the "collec- 
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J = 2 J = 4 J = 6 J = 8 J = 9 

7rd 5/2 1.0 (50%) 2.0 (50%) 2.9 (48%) 3.6 (45%) 3.9 (43%) 

7rd 3/ 2 0.0 (0%) 0.0 (0%) 0-1 (2%) 0.4 (5%) 0.6 (7%) 

vd 5/2 1.0 (50%) 2.0 (50%) 2.9 (48%) 3.6 (45%) 3.9 (43%) 

vd 3/2 0.0 (0%) 0.0 (0%) 0.1 (2%) 0.4 (5%) 0.6 (7%) 

TABLE II: Single-particle components of the total spin in 
the ID cranking calculation. The first two rows correspond to 
proton orbitals, while the last two to neutron orbitals. The 
numbers in brackets are ratios of single-particle spins against 
the total spin. The unit of spin is h. 

205 r 




FIG. 2: Energy surface of the J = 8h state, calculated by 
means of the 3D-cranked HFB method. There are two minima 
as intrinsic states, which are seen at (8,(j>) — (0°,0°) and 
= (90°, 0°). 



tivity" in this nucleus is attributed to gradual excitations 
into the d 3 / 2 orbit. 

Next, let us present the results of the 3D-cranking cal- 
culations. In this paper, we focus on the analysis of the 
J 11 = 8 + state, where its triaxial deformation becomes 
substantial (7 ~ —10°) in the ID-cranking calculation. 
Obviously from the Fig|21 two configurations compete 
with each other. This competition can be considered as 
a kind of "level crossing" between two different states (or 
configurations), but in the 3D-cranking calculation each 
"level" is represented by a curved surface. There are 
mainly two minima in the energy surface for the J = 87i 
state (Fig©: {6,<f>) = (0°,0°) and = (90°, 0°), and they 
characterize the two configurations. The former mini- 
mum corresponds to the ID-cranking solution where tri- 
axiality is calculated to be 7 ~ —10°. The rotation axis 
is found to be along the shortest axis, and the corre- 
sponding state is expected be of low-K character. The 
latter minimum is the energetically lowest state (yrast 
state) at this spin, which is about 2.5 MeV lower than 
the ID-cranking solution. This yrast state at J 71 " = 8 + 
is found to be axially symmetric because the triaxiality 
is calculated to be 7 ~ 0° (/3 ~ 0.24). In this case, 
the rotational axis points along the longest axis of the 



axially symmetric shape, so that the rotation is of single- 
particle character. As a result, the major components of 
the state should be of high- K characters. Studying in de- 
tail the microscopic structure, the total spin is found to 
be constructed almost purely by the d 5 / 2 orbits (in both 
protons and neutrons): 3.97?i each by the d 5 / 2 orbits of 
protons and neutrons. In addition to the difference in the 
deformation, a lack of the d3/ 2 component in the first 8 + 
state implies the yrast state is surely different from rota- 
tional members of the g-band, from a microscopic point 
of view. The shell model calculation by Wiedenh ver, et 
al. 0] says that such a configuration, that is, (d 5 / 2 ) 8 , 
corresponds to the first 8 + state (which does not belong 
to the g-band) observed in experiment. Therefore, our 
result is consistent with the shell model calculation, too. 

From this result, we can conclude that the yrast state 
found in our calculation at J — 8h is a high-it' state 
with K v = 8 + . It was experimentally observed to be 
energetically lower than the second 8 + state by about 2 
MeV. 
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FIG. 3: (Top panel) A cross section of the energy surface at 
= 0°. (Bottom panel) Triaxial deformation at cj> — 0° . 
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0.31 0.25 


0.21 


0.18 


0.14 0.11 


0.17 0.21 


0.23 0.24 



TABLE III: Change of f3 as a function of the tilt angle 9, 
obtained in the 3D-cranking calculation for J = 8h amd <j> = 
0°. 
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In the lower panel of Fig0 the triaxial deformation is 
seen to have 7 ~ 120° at 20° < 6 < 45°. Because the 
"level crossing" happens at 6 ~ 20°, we cannot exactly 
see how the graph continues toward 9 — > 0. However, 
from the trend of the graph, it is possible to guess that the 
graph forms a symmetric shape with respect to 9 = 45° in 
the upper panel of Figg], and that 7 -> 120° for 9^0° 
in the lower panel of FigEl This reflection symmetry 
around 9 = 45° indicates that the state in 0° < 9 < 45° 
and the state in 45° < 9 < 90° have the same intrinsic 
structure. 

CONCLUSION 



that is, through the self-consistent 3D-cranking model. 
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The ground-state rotational band was studied with the 
self-consistent ID-cranking calculation. It was confirmed 
in this study that an effect of triaxiality in the nature of 
the rotational band is important at high spin, as previ- 
ously anticipated by Valor, et al. In addition, two high- 
spin states at J = 8 + observed in experiment were anal- 
ysed by means of the self-consistent 3D-cranking calcu- 
lation. It concludes (in a qualitative manner) that the 
yrast 8 + state is an axially-symmetric high- if state cre- 
ated by the deformation-aligned protons and neutrons in 
the d 5 / 2 orbitals while the second 8 + is a rotational mem- 
ber of the ground-state rotational band with substantial 
triaxial deformation. 

For the first time, in the framework of the self- 
consistent and microscopic method, the two V = 8 + 
states in 24 Mg, which correspond to low- and high-if 
states respectively, are explained on the same footing, 



[1] P. Bonche, H. Flocard, P.-H. Heenen, Nucl. Phys. A 467 
(1987) 115. 

[2] W. Koepf, P. Ring, Phys. Lett. B 212 (1988) 397. 
[3] S. W. Robinson, R.D. Bent, Phys. Rev. 168 (1968) 1266. 
[4] I. Wiedenhover, et al., Phys. Rev. Lett. 87 (2001) 142502. 
[5] R. K. Sheline, I. Ragnarsson, S. Aberg, and A. Watt, J. 

Phys. G 14 (1988) 1201. 
[6] A. Valor, P.-H. Heenen, and P. Bonche, Nucl. Phys. A 

671 (2000) 145. 
[7] M. Oi, P. M. Walker, and A. Ansari, Phys. Lett. B 525 

(2002) 255. 

[8] K. Kumar and M. Baranger, Nucl. Phys. A 110 (1968) 
529; M. Baranger, and K. Kumar, Nucl. Phys. 62 (1965) 
113. 

[9] N. Onishi, Nucl. Phys. A 456 (1986) 279; T. Horibata, 
and N. Onishi, Nucl. Phys. A 596 (1996) 251. 
[10] P. Moller, J. R. Nix, W. D. Meyer, and W. J. Swiatecki, 
At. Data Nucl. Data Tables 59 (1995) 185. 



